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In this work I study the polarization of the Cosmic Microwave Background (CMB) induced by a cosmo- 
^ i- logical background of gravitational waves (here GWs), probably originated during the inflationary epoch. I 

' discuss the influence of a possible reionization of the Universe, which should happen at late times {z ^ 20-;-30): 

I show that, in the presence of such a reionization, we have a remarkable enhancement of the present po- 
' larization of the CMB. I also point out the role of a background of gravitational waves of large wavelength, 

. re-entering the horizon very late, in producing the anisotropy and the polarization of the CMB through the 

Sachs- Wolf effect. Then, I use the standard formalism of the Stokes parameters to describe the properties of 
^ . the CMB photons and I study their evolution with the time, the energy and the spatial coordinates through 

' the radiative Boltzmann transfer equation. I solve such an equation in the form of a second-kind Volterra 

0^ ■ integral equation, by using the analytical method of the iterated kernels: this particular method can be 

successfully applied when we have an exponential kernel, as in the case when the Universe knew a fast and 
. late reionization, due, for instance, to the Hawking evaporation of primordial black holes. 
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1. THE RADIATIVE TRANSFER EQUATION AND THE SACHS- WOLFE EFFECT. 

The equation which expresses the evolution of the Stokes parameters with the time, the energy and the 



$H ' space coordinates is the Boltzmann radiative transfer equation ^ : 
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(rj is the comoving time, defined as J[dt/R{t)], R{t) = being the scale factor of the Universe). 
' Here, the variables are the components of a 4-vector Ua = na{0, </>, v), which is a function of the polar 

angles {9, (f)) and of the photon frequency; these components are the four Stokes parameters which describe 
the polarization status of the radiation ^: 

= {Ii, Ir, U, V). (1.2) 

Ir are the left and right intensities of the radiation, while U and V respectively represent the linear and 
the circular polarization of the CMB photons. In the Chandrasekhar formalism ^ one has: 

I = Il+Ir Q ^ h-Ir. (1.3) 

In eq. (1.1) 7" are the components of an unit vector in the propagation direction of the photons (in 
particular, the z axis is taken parallel to the propagation direction of the gravitational waves we consider), 
(Jt ~ 6.65 X 10^^''' cm? is the Thomson scattering cross-section, R{rj) is the scale factor, /i = cos6', 
P(/i, 0, /i', (/)') is the Chandrasekhar scattering matrix and N,, is the comoving number density of the free 
electrons: such a variable introduces in the Stokes parameters evolution an important dependence on the 
ionization history of the Universe. In fact, the Thomson scattering of an anisotropic radiation (as the CMB 
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is) giv(s it a linear polarization ^: clearly, in a rcionized Universe, the; Thomson scattering processes are 
effective for a longer time, thus one obtains an enhancement of the resulting CMB polarization. 

The relevant term in the left-hand side of eq. (1.1) is the third one; in fact, it expresses a very peculiar 
effect due to a cosmological background of gravitational waves, known as Sachs- Wolfe effect ^ . The presence 
of a ripple in the space-time metric produces a shift in the photons frequencies: in fact, a time-dependent 
perturbation in the metric, i.e. a gravitational wave, determines the non-conservation of the energy along 
the photons sight-line, thus causing a frequency shift: 

dr] 2 dr] ' ' ' ^ ' ' 

If one writes the perturbation of the metric tensor as a plane wave 

ha0 = /I e-'(^--'='')£„/3, (1.5) 

with 

ea/37"7^ = (1 - cos6'2)cos2(/), (1.6) 

the frequency shift becomes 

Thus, the Sachs- Wolfe term obviously influences the evolution of the Stokes parameters and, as a consequence, 
there is a polarization induced in the CMB by cosmological gravitational waves. 

2. THE SECOND-KIND VOLTERRA EQUATION 

Here I will briefly recall the standard formalism of ref. ^, useful to turn the transfer equation (1.1) into 

a second-kind Voltcrra equation. 

By using two new variables, a and /?, and two new vectors: 

a = ^ (1 - m')cos20(1, 1, 0), (2.1) 

b = i((l + ii^)cos2(j), - {I + /i^) cos 20, 4/isin2^!)), (2.2) 
the vector ft is written in the form 

n = no (1, 1, Q) + e~'^^-'^~^'^^ {aa + (3b). (2.3) 

In eq. (2.3) no = (f^Oj^OjO) is the unperturbed solution one has in absence of gravitational waves. 

By expressing /, Q, U as functions of these variables, one obtains two coupled, flrst order differential 
equations ^: 

/? + [q{ri) - ikt4(3= F{r,,k), (2.4) 

i + [q{v) - iMC= H{'n,k), (2.5) 

where q{ri) = aTNeR{rj) depends, through the electron density Nf., on the ionization history of the Universe. 
The right-hand sides of eqs. (2.4) and (2.5) are 



(1 + /x'')'/3(M,fc) - \i{l^,k) (1 - /x'2)2 



(2.6) 
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H{ri,k) = D{k)k 



kr] 



(2.7) 



Eq. (2.6) comes from the calculation of the scattering integral in eq. (1.1), while H{ri, k) is given by the time 
derivative of the metric tensor perturbation, calculated by taking into account the particular gravitational 
waves spectrum predicted by the theory The form of their spectrum is expressed by the factor D{k), 
which is strongly model-dependent: the simplest choice, when quantum gravity effects are not involved, is 



D{k) = Vl67rGi/2(fc); 



(2.8) 



in eq. (2.8) H is the curvature of the de Sitter Universe at the epoch when the GWs enter the horizon. 

More complex forms for D{k), referring to the various kinds of inflationary Universe, can also be con- 
sidered: here, for simplicity, I will discuss the simplest case only. 

By substituting the variables /3 and ^ in eq.(2.6), formally obtained by eqs. (2.4) and (2.5) the 
resulting equation is a 2-nd kind Volterra integral equation which reads: 



F{v,k) = 
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drj F{r] , k) K^(r] 



dri' H{ri\k)e--('^"^ K_{ri ,ri' ,k)y (2.9) 



where 



and 



K±{'n,V ,k) 



"(»?) = / q{v) dri . 



(2.10) 
(2.11) 



is the optical depth. 



3. THE SOLUTION WITH AN EXPONENTIAL KERNEL: 
THE METHOD OF THE ITERATED KERNELS 



In Ref. * I discussed a mixed analytical-numerical method to solve the radiative transfer equation, based 
on the resolvent method ^; such a method can be used when the kernel of the equation has a polynomial 
form. The interested people can find a detailed discussion of this technique in the quoted reference. Here 
I want rather to discuss a totally analytical way to solve the transfer equation, based on the properties 
of the Volterra integral equations. The method of the iterated kernel is a very simple and elegant way to 
solve a Volterra equation, but it can be used in few cases only: the fundamental task is to find a resolvent 
function by writing a convergent series built up with the kernel functions, after they have been convoluted 
and integrated. For a detailed discussion of this method see refs. here, I will briefly recall the main 
results of my analysis, see ref. 

Look now at the general form of a Volterra equation: 

fix) = G{x) + [ K{x,x') fix') dx'. (3.1) 
Jo 

Here K{x,x') is the kernel and G{x) is the source function. Starting from eq. (3.1), the iterated kernel is 
given by ^: 

Km{x,x') = I K{x,t)Km-i{t,x')dt. (3.2) 

If the successive approximations of the kernel obtained in such a way reaUy converge, the integral equation 
can be written in the form 

oo „ 

fix) = G{x) + J2 ^"^ Km{x,x')G{x')dx' . (3.3) 

m=l 
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By defining the resolvent as 

oo 

T{x,x';\) = J2 >^"'~'Krr,{x,x'), (3.4) 

m=l 

the final solution is ^: 

fix) = G{x) + A / G{x')T{x,x':\) dx' . (3.5) 



The applicability of this method clearly requires the convergence of the series in eqs. (3.3) and (3.4); 
moreover, the integral (3.2) should involve quite simple, analytical functions. 

An exponential kernel satisfies these conditions very well: in our specific case, the kernel can be put in 
an exponential form if the function (7(77) (see eq. (2.11)), specifying the ionization history of the Universe, 
has such a form. As I proved in "'12 a 

reionization of the Universe induced by the Hawking evaporation of 
Primordial Black Holes (PBHs) has just an exponential behaviour. 

With this choice, the kernel of the Volterra equation can be approximated by 

S{k,r],r]') = q{r])K+{k,r],r]') ~ exp{ci) exp{c2k{r} - ri )) . (3.6) 

The parameters ci and C2 are obtained by a low fit, whose results are shown in fig. 1 for the case 
z = 20, x = 1, for which ci = 5.22, ca = -18.42. 

Finally, the original Volterra equation (2.9) can be turned 1° into a very simple, linear, first order 
differential equation with a boundary condition F{k, 0) = G{k, 0); the final solution is 

Fik,ri) = G'{k,r,) + - - [ exp{{X + C2k)r,){exp{br,) - 1) J , (3.7) 

where a and b are two constants coming from the calculation of the primitives. 



4. THE CALCULATION OF THE CMB POLARIZATION 



The exact solutions of eqs. (2.4) and (2.5) are obtained by substituting in eq. (2.9) the functions F{r], k) 
and H{rj, k) determined through the analytical solution of the Volterra equation previously found. Then, 
the present CMB polarization and anisotropy, induced by a superposition of fc- fixed GWs, are given by i": 

P\k) = \P{k)f [(1 + + 4^2] d„, (4.1) 

A\k) = £^ \ak)\' (l-n'fdn. (4.2) 
The normalization is carried out by defining the polarization degree as follows: 



Pik) = ^ ^ ; (4.3) 

7' (A:) = |a(fc)|2 {\-i?f + 4/2 + 4/0 \a{k)\{\ - y?f; (4.4) 
here ^ 

^° " exp{hv/kBT) - 1 ^^'^^ 

is the unperturbed intensity; I take also 

hV ~ 10"^^ GeV, ksT ~ 2.35 x 10"^^ GeV. (4.6) 
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In the same way, the normahzed anisotropy is defined as: 

a {k) = -2— , (4.7) 

I (fc) 

2 

where / (fc) is still given by eq. (4.4). 

In fig. 2 I showed the polarization peak obtained from the solution of the transfer equation in the case of 
a standard ionization history of the Universe in figs. 3 and 4 I showed my results respectively in the case 
of a rcionization which happens at z = 10 and at z = 20. I plotted the natural logarithm of the polarization 
degree vs. the wave number of the considered gravitational waves. In fig. 2 the maximum of the polarization 
corresponds to a wave-number k ~ 15, i.e. to gravitational waves re-entering the horizon at the end of the 
recombination epoch: their effect is dominant. The peak is the direct consequence of some mathematical 
suppressions: for rj < 0.03, the source function G{ri, k) is strongly suppressed by the exponential function; 
the physical reason for this suppression is that the polarizing effect cannot be relevant because the Thomson 
scattering has a low efficiency and thus the perturbations arc ovcrally damped. As a consequence, GWs 
re-entering the horizon too early are unimportant in order to produce a CMB polarization. 

For 0.03 < T] < 0.076, Thomson scattering just begins to have some importance, but the polarizing effect 
of GWs having a wavenumber k = 15 is dominant, thus masking that of GWs of different k. Therefore, 
the polarization is relevant only near rj = 0.066. Finally, if none rcionization is considered, for rj > 0.076 
the ionization degree is very small: as a result, the polarizing effect is negligible because we have no more 
free electrons which scatter the CMB photons. 

As one can see from figs. 3 and 4, a sharp secondary peak appears also at smaller /c (A; = 4), in the 
case of a total reionization coming at z = 20. That represents the enhancement of the CMB polarization 
in a reionized Universe: in this case, a polarization for the CMB photons should be observable at angular 
scales near 9°, slightly smaller than the ones tested by COBE. On the contrary, for a standard ionization 
history, one should expect to see the CMB polarization only at angular scales near 2°, thus requesting very 
high precision measurements. 
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